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ABSTRACT 


A layer-ocean model was developed based on the mixed 
layer models of Kraus and Turner (1967) and Denman (1972). 
The ocean model was coupled to Pearson's (1972) time vary- 
ing, symmetrical-stationary hurricane model which was 
based on a model proposed by Riehl (1963), and in which 
the air-sea interaction was specified by Cardone's (1969) 
extension of Blackadar's (1965) saggelionana sy 6X ; baroclinic 
beundary-layer model. Time-dependent solutions for ocean 
mixed layer depth and temperature were obtained in res- 
ponse to interaction with the atmospheric model. Solu- 
tions indicated that the interaction between entrainment 
mixing and upwelling was most important in changing mixed- 
meee depth and temperature. RKadiatidonal effects, inter 
nal waves, turbulent scale energy dissipation, and large- 
scale ocean currents were not included in the present 


ocean model. 
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Lg INTRODUCTION 


ie 1S widely accepted thet the source negqion -foreene 
hurricane heat engine is the underlying ocean. Although 
the ocean is virtually an infinite heat source, it is 
generally accepted that sea-surface temperatures in excess 
of 26C are necessary for forming and maintaining a mature 
hurricane. At the same time the ocean acts as a sink for 
the hurricane's momentum which is ec to the ocean 
through surface stress. It 15s the balance of these two 
processes which was the basis for Riehl's (1963) steady- 
state hurricane model in which the tangential wind profile, 
outside the radius of maximum winds, was specified as 
me! = constant. 

Songqnaty Tia. adeopeca hihehl Ss weriwmeeo ewe lap a 
model which used boundary-layer predictive equations for 
moisture and temperature, and bulk aerodynamic equations 
for the fluxes at the air-sea interface. Corgnati in- 
cluded a fairly simple ocean model which described the 
oceanic Bae KONEC to convective mixing as a result of heat 
extraction; however, most experiments were run with con- 
stant sea surface eas eu eee 

Pearson (1972) modified Corgnati's model by uSing 
Cardone's (1969) extension of Blackadar's (1965) two- 


fayer near=neutral>, baroclinic boundary-layer model)-to 





calculate the fluxes at the air-sea interface. Pearson's 
model produced steady-state solutions for constant sea- 
surface temperatures, and exhibited greater sensitivity 
to these temperatures than did the Corgnati model. 

It was felt that to more realistically describe the 
hurricane-ocean interaction, a more complex ocean model 
incorporating wind mixing, convective mixing, and advec- 
tion, was required. it was the object of this study te 
develop such a model and couple it to Pearson's atmos- 
pheric model through the boundary-layer fluxes of heat 
and momentum. The current ocean model is based on mixed- 
layer models of Kraus and Turner (1967), which ee 
the mixed-layer temperature and depth in response to wind 
mixing and canvection, and Denman (1972), which includes 
the effect of upwelling. It was hoped that the ocean 
model may help to explain which of the mechanisms is of 
greatest importance in producing the areas of low sea- 
surface temperatures observed in the wake of hurricanes. 
Investigators have attributed the low temperatures to 


various mechanisms with little agreement as to which is 


the most important. For example, Jordan (1964) concluded 
that mixing was the more important factor. Fischer (1958) 
considered upwelling to be the most important. Leipper 


(1967) concluded that upwelling, produced by surface 
divergence, was the most important near the hurricane's 
center and that mechanical mixing dominated at larger 


radil. 





The present ocean model was run uSing each of the 
tEheee cooling mechantsmemiCconvect ton, eyand.) maxing wand 
savection), separately and in various combinations, in 
an attempt to study their relative effects upon the 


meeeat Upper thermal strweture and subsequent heat contene. 
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Jf as DESCRIPTION OF STR  MelrL 


joe ATMOSPHERIC (VORTEX 

<a Bo ee oD OP so oe =~ a work, assumed conservation 
Of potential vorticity in the inflow layer and a bulk 
aerodynamic equation for surface Ses caw atlneclonicda aie dia a 


coefficient and inflow angle, to derive the relation 
Vor = constant (1) 


which specifies the radial profile of cyclonic tangential 
winds at the top of the atmospheric boundary layer as 
Sewn in Figure (1)... The outer Te Of thescircu la tao; 
described by Riehl's model was defined by the radius (ro) 
meewhich the cyclonic flow becomes anticyclonic in the 
outflow layer. Riehl assumed that the air ascending in 
the eye-wall region, and flowing cyclonically outward 
aloft, conserves absolute angular momentum. With this 


assumption, one obtains 


(2) 


where f is the Coriolis parameter and the i subscripts 
denote values at the radius of maximum wind speed. The 
present atmospheric model assumed that the radius of 
maximum wind (xr) was coincident with one eye-wall radius. 
Since the eye wall is a preferred region for deep convec- 


tion and latent heat release in observed hurricanes, 


i gal 





those boundary-layer air parcels with the greatest neat 
and moisture content in the model were assumed to be at 
the eye wall. Thus, r, was taken to be the radius at 
which the equivalent potential temperature ere: averaged 
over five grid increments (Ar = 3km), waS a maximum. 
Between the vortex center and r. the winds were assumed 
a 

=u , 

to obey Vox = constant, which agreed with data from 


Riehl (1963), as well as Shea and Gray (1972). 


From observations, Riehl developed the relationship 


P' = -2.56 @ ' (3:) 
€ . 


where P' is the pressure departure from 1005 mb and 8! 


is the equivalent potential temperature departure from 
Beso OK. It was assumed that the atmospheric vortex was 
in a state of quasi-equilibrium with the boundary layer 
fluxes of heat.and moisture, which determined the in- 
crease in 88! and thereby P' of the Cerne aly 
Assuming a gradient balance between the wind and 


pressure field 





2 
1 op 6 
ee eee ve 
Seas . ff Vo (4) 
and using the relation 
Vv 4 ns 
S) i 
a i 

i = = ; C5) 


ah 





integration cf (4) Fon ro te res yields the following 


expression for the maximum wind speed: 


O.5r. 
Vv = —— [(-2fvr.r - Yr, 
0. rours io i 
wey ae Vr Yr Sr) eo ar ee yd (6) 
1 o i fon aul en e 
[ae i O ; 
O 


For the symmetric model, the time-dependent predictive 
equations for potential temperature (8) and specific 
humidity (q), which characterize the boundary-layer pro- 


perties through a specified depth (AP), are 





d(rv_@) gQ 
de, % _ 9(wo) S 
oe Y “ox dP 7 oe + SS te 
a oe peed) Be 5 (8) 
ot r or oP LAP 


where (Q,) and (Q7) are the oceanic sensible and latent 
heat fluxes, respectively. os and L are the specific heat 
and latent heat of vaporization. Radial Cae and vertical 
(Ww) components of motion are related through the con- 


tinuity equation 


ow LL nen (9) 
dP rape) | 
The surface fluxes = oe Ts , ama the integrated radial 
Ss 


motion ie! are determined from the boundary-layer model. 


ies. 








oe BOUNDARY-LAYER MODEL 

Following Pearson (1972), the boundary-layer model 
was a two-layer, baroclinic model developed by Cardone 
(1969) as shown in Figure (2). The Cardone model of the 
marine boundary layer was an extension of the near- 
neutral, fixed-terrain model described by Blacradas (1965) 
The ec used the Vg profile from the atmospheric vortex 
model, together with a computed boundary-layer thermal 
Wind and made use of Monin-Obuhkov similarity theory and 
stationary gee ieee theoyyarro. calculate the £rrction 
Piieerty U,, the stability index L., and the inflow angie 
a. These are then used to calculate a vertical mean radial 
wind, u, WitERine= the surtace ieee eee 
4 


, oF Ze 


aL 
ee oe ee (10) 
1 Oo 


ab 


and the total heat flux (H) from 


3 a 
a + — -_— | 
H = Q. Q., [ Beep a 7 Kan Ci) 
eh: ef2 
where ul= fect ton Velocity, (u7p) , where T is 
the radial stress component, 
- - -4 
26> 6.84 x 10 = +428) ox eh0 : u? —~4743; x 10; 


Uy 


surface roughness (m), 


W(L,)= wind profile function [Cardone (1969)], 
L = Z/L', stability index, 
Z.= heights in the surface layer, 


14 





K = Von Karman Consocant. 


@ = mean potential temperature in the 
boundary layer, and 


a= KL /K petiatbbow of heater ans hier acorns js- 
client to turbulent transfer 
coefficient for momentum. 
The individual On amd o fil aos were Speci redsaccard. 
ing to the Bowen ratio Soe eg) Aare s ere The radial 


wind within the spiral layer was determined through in- 


tegration of the Ekman solution to get 





a ES 37 
= j +—_—_- - 
f Vaz ¥2 Vo sin OSB [cos (a z Bh) 
h 
- sin(at+ =e —J)91) 4) (12) 
where 


(£/2u,Kh)"/7 


ee) 
It 


To include the centripetal acceleration, the Coriolis 

Vo | 
parameter (f) was replaced by (f + o—) The values of 
% in the surface layer (Equation 10) and Ekman layer 
(Equation 12) were then used in determining the vertical 
motion (Equation 9) and the advection of heat and 
moisture [Equations (7) and (8)] in the prediction 


a 


portion of the model. 


CG. OCEAN MODEL 
1. Description 


In an attempt to more realistically describe the 


air-sea interaction associated with a time-varying 


JL: 





hurricane model, it was decided to develop an ocean model 
capable of producing time-varying solutions of mixed- 
layer depth and temperature in response to atmospheric 
meamci ng . Previously the hurricane model produced time- 
Pamving solutions forwconstant sea-suritace temperatures 
until a steady state was reached. 

The present ocean model was initialized with sea- 
surface temperature (Ti), mixed-layer depth ca below 


layer gradient (dT/dZ), and deep-layer temperature (T,) 


d 
as shown in Figure (3). By allowing the ocean thermal 
BimuGcture to change, it was hoped that a more thorough 
understanding of the air-Sea interaction process might 
be possible. The present model is designed for a wind- 
dominated region where there is negligible downward heat 
flux, either by solar radiation or sensible and latent 
heat (condensation) fluxes from the atmosphere to the 
ocean. In the absence of upwelling this precludes any 
posSibility for a decreasing mixed-layer depth. The 
atmospheric forcing of the ocean model was by surface 
Peenetic energy input (wind stress) and upward sensible 
and latent meee (evaporation) fluxes. Both wind Stress 
and heat flux produced changes in mixed-layer depth and 
temperature, surface stress by entrainment mixing and the 
heat fluxes by convective mixing. The entrainment mixing 
was caused by surface generated turbulent motions at the 


aepth of thesmnaxedslayerainterface: 


16 





A partitioning procedure was used to budget the 
total atmospheric stress te) between surface turbulence 
production and generation of a radial mass transport ac- 
cording to Ekman theory. The effects of entrainment 
mixing (turbulence) were assumed 6 be redistributed 
momerormily throughout the layer by turbulent diffusion. 

A radial wind-driven current was Some! from the mass 
meansport and, by employing the continuity equation, areas 
of upwelling and downwelling were calculated. Thus, the 
ocean model incorporated the primary mechanisms believed 
to produce the observed ocean cooling in the wake of 
hurricanes. Internal waves, large-scale currents, and 
turbulent-scale energy dissipation within the mixed layer 
were not included. At six-hour intervals the ocean heat 
content was computed and compared to the initial values 
as a check on the model's solution of mixed-layer depth 
and temperature. 

The time-dependent mixed-layer temperatures 
generated by the Been model produced variations in the 
atmospheric model's Vo Protile, which in turn produced 
Variations in the atmospheric forcing (surface stress and 
heat fluxes). This tended to bring the coupled models 
into a quasi-steady state condition in all cases except 
ones in which advection was included. 

2. Mathematical Development of Ocean Model 
By neglecting molecular heat fluxes and viscous 


generation of heat, Denman (1972) wrote the first law 


17 





of thermodynamics (conservation of thermal energy) for 


an incompressible fluid in the following form: 





Q 
ar J; 
dit: - GORE ae) 
w p 
where 
_ = 4 aoe " 
oF = hee , Specitiies thesheatacource stern, (14) 
p,, = ocean water density, 
= average extinction coefficient, 
2 = specific heat at constant pressure, 
Roeoctdiemadiatwon incident on the sea surface, 
and 
Z = depth in the ocean, positive downward. 
Substitution of (14) into (13) gives the first law in 
ere £Oorm 
aT -YZ 
—=— R 
BIE YRe C15) 
where R = R2/PC 
Pp 
The total derivative was expanded in flux form and (15) 
was written in time-averaged turbulent form: 
= J(u rT) - - ————— 
or, i ve a ral a ae ie (16) 
t 14 ox JZ JZ Y 
or a | sez, d(w'T') 
ac @ VRe er (17) 
where UL = vertical average radial velocity in the mixed 
layer, and T, uu. and w were replaced by Te Ts; ue ai une 
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and, wo+ w' po respectively... The term invediving (ogee) iw 
Equation (16) was dropped because it was assumed that 
there waS no systematic correlation between u' and T'. 
Equation (17) is the heat conservation equation where 
d(w'T')/9Z represents the local divergence of the turbu- 
lent vertical heat flux. It was this term which redis- 
tributed the heat exchanged at the boundaries uniformly 
throughout the homogeneous mixed layer. At the ocean 
surface, the downward turbulent heat flux, (CE was 
equal to the nonradiative heat transfer through the 


ocean surface, i.e., 
(w'T i = - (Q,, + OF Ao (18) 


where oe and Q. were the upward latent and sensible heat 
fluxes, respectively, computed in the boundary-layer 
model. At the bottom of the mixed layer, the downward 
turbulent heat flux was equal to the heat transfer through 


the interface due to entrainment mixing, i.e., 





dh. = 
w ! St = =a tS ee be i 
( dn a ) ( S h? Oe), 
where dh | 
rrr ea time variation of mixed-layer depth, 
w = vertical advection velocity at depth (hod 
Ts = mixed-layer temperature, and 
T, = temperature at a finite depth below mixed- 
h 
layer depth. 
Hor Entrainment mixing, the turbulent flux (wit), was 
O 


positive which carresponds to warm water (T'> 0) being 


ibe, 





forced downward (w'> 0) and cold water (T’< 0) being 
lifted upward (w'< 0). This process was accomplished 
against buoyancy forces by the kinetic energy input at 
ne surface and the potential to kinetic energy conver- 
sion by convective overturning. 

Integrating, Equation~({17) over the depth of the 
Mixed layer (ho) and replacing the boundary turbulent 
heat fluxes by Equations (18) and (19), yields a thermal 
energy equation for the mixed layer, 


aT 


a) x sD cs _ (ah _ 5 
hy = R(d-e 1°) -(0,400)/0,C,- (Ge -w) (Pet 


ae ) (20) 


rf 
By expanding the total derivatives in flux form and ig- 
Metmrng the radiation effects, #R), (20) becomes a predic- 
tive equation for the local change in the mixed-layer 
depth (ho), Le s:; 
dh ~h oT a (ier ) d(w T_) 
jie 1 Ce ee 2 es ee ees 
dt Chea 











+ ow. (2) 


To develop a corresponding equation for mixed-layer tem- 
perature changes, Denman (1972) integrated Equation (17) 
twice, once over a depth z and once over the depth of the 


mixed layer (ho). The resulting expression, 


20 








Qu 
II 





az (22) 


Qu 
ct 


represents a mechanical energy equation for the mixed 
h ee 
ilaver. The term JS a (w'T')dZ specifies the conversion 
O 
Siepotential energy into kinetic energy by convection 


within the mixed layer. If one considers the mechanical 


energy balance expression 


Wot G =) D.= 0 (23) 


where G = kinetic input from the wind, 
D = dissipation within the mixed layer, and- 
ho ee 
w= f iv a2 (potential) tomeinetic emery 
e) Convers Lon) >| 


substitution of Equation (22) into (23) yields 


dT oh (0O_+ 0.) 
Se) Se eae 2 . 
aE 7 eC ho hea my aae (24) 


w P 


Expanding the total derivative of temperature in flux 


form and solving for the local temperature change yields 


¢ 


ce. _ 2(g-p) — 7 '2@g7 Qs) 
dt 2 h 
h. . ane 
ia, d(u rT.) ; 0 (w To) se 
r or OZ 


which is a mechanical energy equation specifying the 


change in mixed-layer temperature due to entrainment 


Za 





musing, Convection, and advection. BGUattone( 29) Can 


be used to eliminate oT /dt from the heat equation (21): 


oh ees) (Q,4 Q,) ate 


= ee = -Y 
) a Debts) r or 


+ w (26) 








Equations (25) and (26) represent a system of two equa- 
Erons with three unknowns, T , noe and Th (u and w will 
be obtained from Ekman theory and the equation Of seone 
m@rnaggaty). To close this system therefore, a third equa- 
tion or relation was required which specified changes in 
Th as being dependent upon changes in no: The parameter 
which imposed this dependency was the below-layer gradient 
(dT/dOZ) which was assumed to be uniformly affected by ver- 


mmo! Velocities and thus remained constant. Time varia- 


tions of Th were specified by the expression 


oT dh 





h = Oe. OF 
— = - - ~— ad 
at w= at | Beh . ae 
O 
where w = vertical advection Ve Woerty . 


The solution for T which represented the lower tempera- 


h! 
eare of Mer eneratnea water, proved to be very eruciatl 

in the mixed-layer depth and temperature solutions. at 
was the temperature difference (T= TY which imposed the 
primary retarding force against unlimited deepening. As 


time progressed the difference (T z increased, forming 


a 


-teeacetal Step in the temperature profile at the base of 


ZZ 





the mixed layer. Because of the stable lapse-rate 
(d9T/OZ < O), upwelling (w<0O) contributed to a lowering 
of Th and, hence,».sincreased thesditftienrence es Th): 
Thus, upwelling had the effect of retarding the tendency 
for mixed-layer deepening. At the same time upwelling 
(w<O0) enhanced mixed~layer temperature changes by de- 
Sereasing the hayes: thickness (Equation 26) so that the 
energy input and heat fluxes would have a greater effect 
Ogeuatlion 25 )u Thucem Bawatrons (2598826) , and (27) 
formed a system for unknowns ee Toe and Th Space and 
time derivatives in (25) and (26) were approximated by 
centered differences, except for the initial time step 
which used a forward difference. A forward step was 
also used every 17th time step (3 hours) to prevent di- 
verging solutions. Forward time differences were used 
in (27). The difference forms of (26) and (27) were 


solved together for he and T. using a £Lirst guess depth 


h 
ie tempera tum e (tT) from the previous time step. The 
first guess depth in (26) and (27) wasS increased by 
0.05m increments Saal the difference between the right 
and left hand sides of the difference form of (26) was 
less than or equal to O.lm. Using the new mixed-layer 
depth, mixed-layer temperature was calculated from (25). 
The mechanical energy input (G), in the ocean 


predictive equations (25) and (26), was specified 


according to Kraus and Turner (1965) as 


ZnS 





p 


I a. 57 2 
G = egal ae / ls (28) 
g os aoe 
where Gg “}=""qravitational awcecleration, 
dp 
il W ad 
C= 2-— ——_ tt, COC erent Of “expansion, 
OL acer 
Ww 
Pp. = alr density, 
oe = sea water density, 
Ue Sse nospher eet r tet ronevelocityjacand 
tA = atmospheric surface stress. 


In the ocean model, surface stress (TO) used for the 
energy input term (G) depended exponentially upon the 
ratio of the mixed-layer denth (ho) and the roughness 


length (20): 


t= Lt. EXP =-2. (29) 
a 


where 


T = total stress available from boundary-layer 
5 model, and 
C = an arbitrary constant. 
h 
For a given roughness length (Zo), the function exp Gee 
| O 


specified the percentage of total stress (T_): available 
for turbulence generation as a function of mixed-layer 
depth (Table I). The use of the above partitioning pro- 
cedure was an attempt to budget the available energy 
between turbulence generation (G) and production of wind- 


driven currents in the ocean model. Experiments were run 
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with various values of C and compared Eo veachmotinc. ea. 
well as to a run in which all the surface stress was 
Mellized for turbulencewgeneration. The dissipation .p: 
was set to zero in all aie nce 

The radial velocity, ue which is needed to close 
the system, was calculated in accordance with Ekman theory 
which specifies a spiral layer to a depth (2) where the 
current direction has reversed from its surface direction 
and the speed is exp (-7) of the surface speed. Below 
the Ekman layer it was assumed that the wind had no direct 
ioe limence . ACGOraing to Ekimanmethe radial component Jon 


motion within the spiral layer was described by 


(T -T )cos a 
ann © 


uu, = oo ee EXP (-BZ) [cos (~BZ)- sin(-BZ) } (3:08) 

pi (2e oe 

W 

where Z = 1/B, specifies the Ekman depth, 
i ly 2 

B = || ’ 

2U/) : 
f = coriolis parameter, 
om = sea water density, 
Hw = vertical eddy coefficient of viscosity, 
Z = depth, and 
a = atmospheric inflow angle. 


Integration of (30) over a specified Ekman depth (2) 
yields an expression for the net mass transport in the 


madial direction for the Ekman layer, 


a = ny (31) 
r 8 
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where Tp = (as T.)cos QQ; tangentlalvyatmospheric. siumeace 
stress. ' 


This represents a radial mass transport per unit distance 
mm ene tangential directron. One can calculate a mean 
Padial advection velocity (a) For the mixed layer from 


the relation 
u Z = M7), : (32) 


The mean radial (u_) and vertical (w) COmponents OlLemot.on 


Mare related through the continuity equation 
d(p SO 6) 


d(p_ w) 
Wo ee (33) 


OZ Yr or 


Equation (33) was integrated over depth (2) using the 





boundary condition w(z=0) = 0, to obtain wat Z 
/ boy 8 
= di, 0 0 
ee rs 34 
p., W Soa ee ) (34) 
which was taken to be w at nos The probable occurrence of 


pemeecurn flow at some depth below the mixed layer was not 
included in the present ocean model. Thus, the only effect 
on the Ipeienelasreiz gradient (dT/dZ) was that due to verti- 
cal velocities and, hence, the gradient remained constant. 
As a check on Eivelic ocibina processes (wind mixing, 
convective mixing, and advection), the initial ocean heat 
Gontent (epths,/ crime for a depth Z, = 130m was calculated 


b 


mesed on the initial temperature profile (Figure 3). Inxs 


a 





value was then compared to the heat content at ecachecr ima 
peint, computed at G—-hour intervals. ~FPoOr=econvecctive wr: = 
ing-only, with the kinetic energy input (G) and the 
advection velocities (a, and a icoccel te Zero, Hoavatiens 
(25) and (26) predicted excessive mixed-layer deepening 
mmaecoOling. Also the final grid point heat contents did 
mor reflect just the heat (H) lost to the atmosphere due 
to upward heat fluxes coe and i Hence, for this case 
the mixed-layer depth and temperature were specified by a 
miemmal equilibrium problem using the following set of 


equations: 


Pie Te. 
s os 
= t | ees "y+ 
Eee Ts UO A ogc an. Am 
h = h ' + Ah (35) 
O O 
oT 
T = T ' + = Ah 
Ss Ss 0Z 
where hove initial mixed-laycigdepth, 
ee initial mixed-layer temperature, 
oT/OZ. = below-layer gradient, 


oe 2.7 Os POedike Wey aiscee me ae oc, 
p= sea water density, 

At = time step (sec), and 

= specific heat 


For convective mixing only, equations (35) were solved at 


each time step for the change in mixed-layer depth (Ah). 


With this calculated depth change, the new layer depth (ho. 


and temperature (tT) were specified. 
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Figure 1. Atmospheric vortex model. 
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Figure 2. Boundary layer model. 
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ile DISCUSS ON WOOF estes 


As stated previously, the Siyect ew thas os bucdy™= was tO 
develop an ocean model capable of producing time-dependent 
solutions for mixed-layer depth and temperature in response 
to forcing from a symmetrical, stationary hurricane model. 
Previous work with the hurricane model was done with con- 
stant sea-Surface temperatures and, consequently, little 
was learned of the ocean thermal changes or of the corres- 
ponding changes produced in the hurricane model. 

In this investigation three mechanisms believed to con- 
tribute most in changing the ocean thermal structure were 
isolated and studied. These were entrainment mixing across 
a stable layer (wind stress), convective mixing (latent 
and sensible heat fluxes), and vertical and horizontal 
advection. It was hoped that the combined time-dependent, 
atmosphere-ocean model would help determine which of these 
mechanisms, cee acting alone or in combination, was most 
important in preducing the cold sea-surface temperatures 
observed penegtn hurricanes. To accomplish this, a series 
Of experiments were run uSing each of the cooling mechanisms 
alone and in combination with each other. Figures (4) 
werough (32) depict radial profiles of various ocean para- 
meters for the different experiments. For these runs the 


ocean model was initialized with mixed-layer depth of 


on 





30 meters, mixed-layer temperature of 30C, below-layer 
Gradient of 10C/100m, andw@deep=ayer Wemperature of BOC. 
The lower boundary of the model was fixed at 130 meters 
(Figure 3), below which there were assumed to be no 
thermal effects. Other initial temperature profiles can 
be used; however, the above values were believed to be 
fairly representative of the tropical ocean regions. The 
atmospheric model was initialized with air temperature of 
29C, radius of maximum wind (xr. = 27km), and maximum wind 


(V5 = 30m/sec). In runs with various constant sea-surface 


af 
temperatures, Pearson (1972) found a steady-state v Of 


a 

60m/sec for 30C. 

Mie vaiue of the constant (C) in ene Wind stress parti- 
tioning equation (29) was arbitrary fixed at no in all 
the initial runs except the advection-only case. For this 
run the partitioning process was not used; consequently, 
all the available stress ces) was used to produce radial 
Wind-driven currents (u_) and vertical velocitwres (w). 
Since the nae tient nt process was an important assumption 
in the model, an entrainment-convection-advection experi- 
ment was — With a value of 5 x one for thervcons tamtegicc) 
imenauation (29). With a smaller valwe of "C; Equation (29) 
partitioned more atmospheric stress to producing ocean 
currents and less to entrainment mixing at all mixed-layer 


depths (Table I). An entrainment-convection run was also 


made in which the partitioning process was not used, hence 


4 





all available stress from the boundary-layer model was 
used in the kinetic energy term (G) in the mixed-layer 


gepth and temperature predictive Equations “[((25) and (2G 7 


Ds CONVECTIVE MIXING ONLY 

For the convective mixing—-only Case, Equations (35) 
were solved for the mixed-layer depth (eee 4) and tem- 
perature (Figure 5). The atmospheric forcing was by 
latent and sensible heat fluxes from the ocean surface to 
the atmosphere as computed in the boundary-layer model. 
Maximum cooling (0.5C) and deepening (5 meters) of the 
meecea layer at 18 hours occurred under the region of maxi- 
mum wind (r.) as expected from Equation (11) which relates 
Octal oceanic heat loss (H) to the cube of the friction 
fmeecity (u,). Oceanic heat loss computed from initial 
and final temperature-depth profiles (Figure 6) in the 
ocean model compared guite favorably with the total heat 
extraction (H) computed in the boundary-layer model 
Ceagure 7). PiewaAtoOuntsOt OCeanwecoolinguduesto ‘eonvece lon 


was comparable to Jordan's (1964) estimates. 


B. ENTRAINMENT MIXING 
Radial profiles of mixed-layer depth and temperature 


due only to entrainment mixing (Figures 8 and 9) were ob- 


moaned by iterative solutions of Equations (25), (26), 
and (27). Thememwas a0 c@ean cunmeniPeqencratie nner 
hun, and the partitioned stressmz9), with C = Ton was 
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the forcing which produced entrainment mixing through sur- 
Pace turbulent generation (G)- Maximum ‘Coo lineml0. oC) 

and deepening (23 meters) at six hours and 1.5C and 45 
meters at 18 hours again occurred beneath the region of 
maximum kinetic energy input (28) as expected from (25) 
ama (26) . The increased cooling and deepening with time 
at larger radial distances was due to the outward migra- 
tion of the hurricane eye wall (r.) in response te vcolder 
sea-surface temperatures near the vortex center. 

Since entrainment mixing produced see vertical heat 
fmoeeaeriigntiion, the heat content of the ocean beneath the 
atmospheric vortex should be conserved except for numeri- 
wore errors. rreure. (10) depots the radial distriscutioen 
of oceanic heat loss at 18 hours for the entrainment-only 
experiment. The area averaged heat loss for the solution 
Gemieain (Table II) was 137 orev enae Considering the 
sensitivity of heat content to changes in temperature and 
depth and since the depth and temperature solutions were 
computed numerically, and are subject to truncation errors, 
mime values in Figure (10) are considered within acceptable 
limits. For example, in the region of maximum mixed-layer 
depth (75m), a temperature difference of 0.1C will give 
a heat contant change of 750 eae The maximum error 
in Figure (10) is 800 SAS - 

Figure (ll) depicts ocean temperature profiles at 6, 


2, and lLé hours as computed in the model. Note the strong 
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Figure 10. 
ocean heat loss for entrainment mixing. 
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gradient which develops below the mixed layer. As men- 
ErtOoned in a previous section, the temperature diftterenece 


ie tT) is the primary retarding force whienwbalancec 


oT 
OZ 


iemedh important factor in determining the depth at which 


piiematmospheric forcing: The below-layer gradient (-~—) 
ease balance will occur. FOr “a (91 Ven energy pu emcee 
one expects deeper mixing for a weak gradient and 


shallower mixing for a strong gradient. 


Cus ADVECTION-ONLY CASE 

For the advection-only case (G=0O, and 27 Nee: mixed 
layer depth and temperature changes were determined by the 
geavection terms in the predictive equations [(25) and (26)]. 
iirempartitioning process (29) was not used, hence, all 
available atmospheric stress re was used to produce a 
radial mass transport (32) and a mean radial wind-driven 
Sumerent (33). The rapid increase in UL (Figure 12) from 
the vortex center to the radius of maximum er (x) pao 
duced a region Of Nhowizontal divergence “and, hence, enor] 
IMe@usly strong upwelling (Figure 13). Just beyond radius 
rs AP emali wegron of convergence and, hence, downwelling 
occurred. However, the gradual decrease in UL was 
balanced by the increased radius EO Matlmtain alrate ly Seon 
Seahe COtal mass transport. Consequently, areas of signi- 
ficant downwelling were not evident within one solution 
Gomain (300°>°km) . 

The oP ROGES of strong upwelling were quite evident in 


Phewmtuxcd— layer Geptn and temperature proriles (Pigurce 
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14 and 15). Within six hours the Maxed Slaver, eneaeeen 
vortex center, was eliminated and the sea-surface tempera- 
ture decreased until the cold deep ocean water (20C) was 
wpwelled to the surface. At 18 hours the radius of maxi- 
mum wind (r.) had moved outward in response to cold sea- 
Suctace temperatures at the vortex center (Figure 16). 

The region of maximum upwelling likewise moved outward, 
leaving a broad region of 20C sea-surface temperatures 
between 15 and 45 km and zero mixed-layer depth to a 
radius of 95 km. The dramatic effects of advection were 
also evidenced by significant ocean heat content changes 
(Figure 17). The positive values correspond to heat loss 
regions and the negative values to heat gain regions. 
Vertical and radial velocities produced a redistribution 
of heat by changing the ocean temperature profile (Figure 
18). Because of the ocean's response to the atmospheric 
stress profile, radial currents (a) transported heat 
across the model's open boundary at 300 km. Consequently, 
there's a significant averaged heat loss (Table II) in 


the ocean model for this experiment. 


DD. ENTRAINMENT AND CONVECTIVE MIXING 

This experiment was similar to the entrainment-only 
mum except heat £lux (Q + 2) was included in (25) and 
ico) os For the initial run, Stress Ieper ettreonmitng (2.7) eva 
utilized with C = Woes Radial profiles of mixed-layer 


depth (Pigure 19) and temperature (Piguve 20)metosely 
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resemble the corresponding profiles for entrainment—only 
(Pagures 8 and 9). The inclusion cf the HReatwE las secan 
produced a slight increase in the amount of deepening and 
eooling compared to entrainment-only; however, the amount 
meas less than 2 meters and 0.5C at 18 hours. One may 
conclude that entrainment mixing dominates when compared 
Memconvective mixing for the initial conditions and atmos— 
fameric forcing utilized in the model. 

Bor this run the oceanic heat content should net be 
conserved because of the latent and ean eros heat exchange 
with the atmosphere. Figure (21) depicts radial profiles 
Semoceanic heat loss, calculated in the ocean model, and 
heat extraction (H) from the boundary-layer model. The 
difference between the profiles represents the loss attri- 
pieced tO inaccuracies in the numerical solution of (25), 
mo), and (27), and Re or une ae 1s as large as H itself. 
Mme vertical ocean temperature profiles (Figure 22) were 
eamrlar to those for entrainment-only (Figure 11). For 
comparison, Figures (23) and (24) show radiai protiles "ef 
mixed-layer depth and temperature, respectively, for an 
Mmeeeimenit convection run in which stress partitioning 
was not used. All available stress oe from the boundary- 
layer model was used in ene Kinetle energy term —(G) tet 
(25) and (26) to produce changes in mixed-layer depth and 
temperature. The increased atmospheric forcing was 


evidenced by an increased maximum mixed-layexr depth (12 m 
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Qeeper) and a lower minimum mixed-layer temperature (.5C 
Golaer) thangror the partrter1oneda run. In both of “these 
entrainment-convection runs, as well as the entrainment- 
only run, the amount of sea-Surface cooling near the vor- 
tex center was sufficient to cause the radius of maximum 
wind (x) to migrate outward. The intensification of the 
atmospheric forcing at larger radii produced a correspond- 
ing increased deepening and cooling of the mixed layer at 
larger radii. Since there was no upwelling in these runs, 
the radial profiles of mixed-layer depth and temperature 
were not modified, and the radius of maximum cooling and 
deepening remained at 20 km. 

In the next set of experiments eeisn will be in- 
cluded and the effect of variations in the partitioning 


constant (C) will be examined. 


EK. ENTRAINMENT, CONVECTION, AND ADVECTION 

This run incorporated all the mechanisms believed to 
be of importance in cooling the ocean beneath a hurricane. 
The total available atmospheric stress (tT) was partitioned 
m9) into that which produced surface turbulence (To) and 
that which generated wind-driven currents Wa EOE 
mae initial run the partitioning constant (C) in Equation 
(29) was set equal to Tes. This gave a stress vs. depth 
profile as depicted in Table I. In addition to ,theskinerrc 
energy and flux terms, advective terms were included in the 


SOlution fou mixed—-layer depth (26) and temperature (25). 


STS 






The radial advection velocity (a) was Specified according 
=O (38) and areas of vertical velocities (w) were calculated 
feom (34). Radial profiles of u (Figure 25) and w (Figure 
26) were similar in nature to those for advection only 
Wagures 12 and 13). However, Since stress partitioning was 
not used for the advection-only experiment, the magnitudes 


of uy and w were smaller in the current experiment. 
As in previous runs, the hurricane eye wall este: out-— 
ward in response to cold sea-surface temperatures near the 
memtex center (Figure 16). The GEueniane region (Fagure 26) 
likewise moved outward in response to the movement of maxi- 
mum surface stress. The changes in mixed-layer depth and 
temperature produced by entrainment and convective Mixing 
(Figures 19 and 20) were significantly modified by the out- 
ward migration of the upwelled region. For the entrainment- 
convection mixing case the region of maximum mixed-layer 
deepening was coincident with the region of maximum cooling. 
However, with the addition of advection, the maximum deepen- 
mee (Figure 27) Syomeneeee! at radii just outside the radius 
of maximum wind (r.), where upwelling was weak and the atmos- 
pheric Porenis (G and Q..+ Q,) was fairly large. At smaller 
radii the upwelling dominated the effect of mixing, leaving 
very shallow Maes Naver Gescne: The region of maximum 
ocean cooling (Figure 28), which remained in the vicinity 
of 15 km, represented an interaction between mixing and ad- 


Wect LON. By decreasing the layer thickness, upwelling 
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enhanced the effects of atmospheric forcing on the mixed- 
ihayer “coolLung) (25) Stress: partitioning iInereasedweme 
percentage of total stress used in the kinetic energy term 
(G) as the mixed-layer depth decreased. There was a corres- 
ponding decrease in the stress used for currents (upwelling) 
which tended to keep the mixed-layer depth from becoming 
zero. Since ultimately it was the value of the constant 
(C) which determined whether mixing or upwelling would have 
the greatest influence on the mixed-layer depth, the entire 
Moermulation of stress partitioning deserves considerable 
Paewe study. 

Smeaniticant ocean heat loss (Figure 29) occurred inside 
30 km Rivets upwelling had the greatest effect on the verti- 
fer temperature profile (Figure 30). At larger radii, where 
entrainment and convection dominated, ocean heat content 
changes reflected the heat extraction (H) by the boundary- 
layer model. wis area averaged heat loss (Table II) for 
the solution domain was 2300 ca leyene at FS" hours. “ AS" in 
the advection-only case, this loss reflects the transport 
of heat across the model's open boundary at 300 km. 

me test the influence of Stréssmpartwtioning (29), a 
second experiment was run uSing a value of C = 5 x one 
With this smaller value, the percentage of total stress 
a) used for entrainment mixing was smaller and decreased 
more rapidly with mixed-layer depth and, consequently, the 
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Figure 30. Same as Figure 11 except entrainment, 


convection, and advection. 


5)2, 


was larger and increased more rapidly with mixed-layer 
depth (Table I). Comparing radial profiles of mixed= 

layer depth and temperature for C = 10° (Figures (27) and 
3) | Swath those for C = 5 x 10° [Figures (31) “amd (32)), 
memes apparent that advection by horizontal sand vertical 
wmelOocities had a greater influence in the latter case. 
[Note similarity between Figures (31) and (32) and those 
Mem mdvyection only Figures (14) and (15)]. For € = 5 x 10° 
upwelling has eliminated the mixed layer in a small region 


near the vortex center and the combined effects of upwell- 


ing and mixing have cooled the sea-surface temperature to 


20C. At larger radii where mixing dominates, there was 
, a 4 
G@ieater deepening and cooling with C = 10 than with 
3 


@e- 5 x 10 . 

Table II summarizes some significant atmospheric and 
ocean variables. The air-sea interaction was evident in 
all cases where there was pronounced cooling of the ocean 
surface near the vortex center. The lower sea-surface 
temperatures reduced the air-sea temperature difference, 
Mence the heat fluxes from the ocean to the boundary layer 
were reduced. The tendency for the model's eye wall to 
seek a region of maximum equivalent potential temperature 
(8) caused the radius of maximum wind (x5) EO Move CWUEWalea 
as the inner surface was cooled (Figure 16). This outward 
migration maintained a sufficient equivalent potential 


temperature difference cr - om ) and, consequently, the 
rk e) 
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Pargure 31. Same as Figure 27 except C=5x10 . 
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maximum wind (Vp ) remained fairly constant. Sebecausewo: 
Nye 7 | 

the Var = constant iprotile, the outward m1 Geate ome. 

r. produced an unréealirstic Mntensitication of “the feimcmra— 

fron at large radii. This in turn maintained a constant 

total ocean mass transport with radius which resulted in 

meat icant amounts of heat being transported out the 

model's 300 km boundary to a region where downwelling would 

presumably have occurred. In order to preperly model the 

2 
movement of the eye wall the Vag tt = conStant profile must 


be modified to prevent this unrealistic intensification at 


large radii. 
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TABLE 2 


SHgOnILficant atmosphere and oceamc 
variables for different experiments 
at TS sheuwies:. 
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TEMES SUMMARY AND CONCLUSIONS 


Mvewresulvts of sthacsvstudy indweate the feasmiarlity “of 
eoupling an ocean modeleto an atmospheric vertex model to 
simulate time-dependent air-sea interaction. Since the 
mee] wasestatdwenary, andwthe effeets of large-scale oe@ean 
current, turbulent-scale energy dissipation, radiation, 
and internal waves were not included, the magnitude of the 
ocean thermal changes may not be representative. However, 
the relative changes for the various cooling mechanisms 
(convective mixing, entrainment mixing, and advection) were 
meparent. 

BU@st, GOnvective miximg hadsdinsigoniisveant effects on 
an ocean thermal structure initialized with 30 meter mixed- 
layer depth, 30C mixed-layer temperature, and below-layer 
gradient of 10C/100m. Convective cooling would have been 
more pronounced for shallower layer depths and stronger 
below-layer gradients; however, these were not tested in 
the model. 

Second, entrainment (wind) mixing produced fairly "subp- 
stantial mixed-layer deepening and cooling in the region 
of maximum wind stress. The radius of maximum wind res- 
ponded to the colder sea-surface temperatures by migrating 
outward. The effect of entrainment mixing, likewise, 
moved outward producing deeper layer depths and colder Sea-— 


surface temperatures at larger radii. The depth and 
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temperature changes due to entrainment and convectleimver. 
comparable to those for entrainment only indicating en- 
trainment to be the dominate process. 

Themed; the effecis of Sere wintdeami ven  cimeulatton, 
when acting alone, produced the most pronounced changes 
in mixed-layer depth and temperature. - Upwelling inside 
the radius of maximum wind du spilaced the relatively warm 
mixed-layer water by cold deep-ocean water advected from 
below. The warm surface-layer water was advected outward 
Py) radiawm cwrrents. Since no significant downwelling 
S@@eurred within the solution domain, heat was lost through 
the model's outer boundary. 

In reality the various mechanisms studied do not act 
separately but rather interact in some complex fashion to 
produce the thermal changes which observations have shown. 
It appeared to be a guestion between entrainment mixing 
and advection as to which has the dominate role in cooling 
sea-surface temperatures beneath a hurricane. Advection 
@efinitely produced significant cooling inside the radius 
of maximum wind but to think of a wind-driven circulation 
with no eae turbulence seems unrealistic. It appeared 
that the interaction of the two processes, with a budget- 
ing of the total atmosphemme stress between turbulence 
generation and current generation, was the dominate 


mechanism for sea-surface cooling. 
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